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Pion condensation in the two— flavor chiral quark model at finite baryochemical 

potential 
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Pion condensation is studied at one-loop level and nonzero baryochemical potential in the frame- 
work of two flavor constituent quark model using the one-loop level optimized perturbation theory 
for the resummation of the perturbative series. A Landau type of analysis is presented for the inves- 
tigation of the phase boundary between the pion condensed/non-condensed phases. The statement 
that the condensation starts at fj,i = is slightly modified by one-loop corrections. The second 
order critical surface is determined and analysed in the /^i — /ib — T space. The dependence of 
the one-loop level charged pion pole masses is also studied. 
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INTRODUCTION 



From the early 70 's, when the idea of pion condensation was first suggested many investigations have been 
performed in the field. Right in the early works, the possibility emerged that pion condensation might have a 
I significance in the evolution of compact stars Beside compact stars, pion condensation can occur in asymmetric 

II ■ nuclear matter as well as in heavy ion collisions at intermediate energies. While the phases of QCD can be well 
Qh' described by perturbation theory at extreme high temperatures and densities, at moderate temperatures and densities 
^ , a perturbative description is not possible, here other methods like lattice field theory and effective field theories can 



be applied. In lattice field theory a second order phase transition was found from normal phase to pion condensed 



se to pi 

phase in two fiavor QCD at a critical isospin chemical potential, which is about the pion mass [1, @, [3|. Pion and in 
^-H ■ some cases kaon condensation was also investigated in numerous effective models such as chiral perturbation th eory 
[8, 'q'I, ladder QCD [lo| . random matrix method ^^l; NJL model (in the mean field approximation) (l^. [T3. [T3. fl5|. 
PNJL model and linear sigma model These investigations mostly focus on properties of asymmetric 

I hadronic matter, however if one wants to describe for instance a neutron star, which is electrically neutral on av erag e, 
. then the condition of charge neutrality as well as /3-equilibrium must be imposed explicitly as discussed in [3, l2Cll |. 
Our main goal here is to go beyond the present status of the field and determine the phase boundary between 
the pion condensed and non-condensed phases in the m — fiB — T space at one-loop level and to calculate the in 
dependence of the charged pion masses beyond tree level. This work is the continuation of our previous works [2ll.[22|. 
[ in the sense that we use the same framework of constituent quark model and the optimized perturbation theory (QPT) 
[2^ 1 as a resummation technique, which preserves the perturbative renormalizability as well as the symmetries of the 
model. However there we used three fiavor model, while the present investigation deals with an effective model of 
2-fiavor QCD. This approximation is sufficient to see the main properties of pion condensation. 
' The paper is organized as follows. In Section II, we introduce the model and discuss its renormalization in detail in 
5^ , case when not only scalar but also pseudoscalar (pion) condensate is present. We treat bosons at one-loop level, while 
the constituent quarks are considered at tree level. In Section III, we present the one-loop equations, which describe 
our system, namely the equations of state for the scalar and pion condensates and the equation for the resummed 
pion mass. Moreover we show that a Ward identity holds, from which the Goldstone theorem also follows. Then we 
expand our equations in powers of the pion condensate, with which we restrict ourselves to a Landau-type analysis of 
the phase transition. At the end of this section a simple formula for the boundary of the pion condensation domain 
is established. In Section IV we diagonalize the mixed boson and fermion propagators, which are determined also 
for small values of the pion condensate. Section V is dedicated to the parameterization of the model, i.e. we fix the 
parameters of the Lagrangian at zero temperature and vanishing chemical potentials at one-loop level. Here we also 
discuss the choice of the renormalization scales. In Section VI we present the numerical results and finally we conclude 
in Section VII. 
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II. THE MODEL AND ITS RENORMALIZATION 



Our starting point is the renormalized SU{2)l x SU{2)ii symmetric Lagrangian with expHcit symmetry breaking 
term, which contains meson as well as quark fields [131 ■ Without the symmetry breaking term the symmetry group 
of the Hght meson sector (tt, a) is 0(4), and the Lagrangian is 



£ = i 



A 



(1) 



where ■0 = ("i dy are the SU{2) doublet quark fields, 4> — (0o, 0i, 02, ^s) = (f, tti, 7r2, tts) are the scalar sigma 
and pseudoscalar pion fields, h is the symmetry breaking external field, and Ta = {tqjITj^), a = . . .3 are the 
quark-boson coupHng matrices. In the Lagrangian above the couplings m^, A and the field are renormaHzed (finite) 
quantities, while 5m?, 5\, and 5Z are the corresponding (infinite) counterterms. Since we treat quarks only at tree 
level no wave function renormalization constant is necessary for the fermionic fields. 

We use the grand canonical generating functional of the renormalized Lagrangian in order to introduce the baryon 
and isospin chemical potentials as well as the scalar and pion condensates. As one shall see, one consequence of the 
presence of pion condensation amounts to the appearance of the non-trivial wave function renormalization constant in 
the equation of state of the condensate. At finite temperature and densities the grand canonical generating functional, 
which determines the A^-point functions of the fields, is 



P0X'ni?-0I?0 Gxp 



where the Hamiltonian of the system is 



dt I d^a;(n0 + 1070^- - W + MbQb + wQi) 



(2) 



n 



1 



^ (tf + (V0)2 + m202) + ^04 _ /i0o + 1^^7,9,^ + ^-^^TM 



and the canonical momenta of the scalar fields are defined by 

6C 



{i = 1,2,3) 



n 



(1 + (5Z)0. 



(3) 



(4) 



In ^ Qb, Qi are the conserved baryon and isospin charges defined by the zeroth component of the corresponding 
Noether currents of the symmetric Lagrangian 



Qb = j A^x^{u^u + d'^d), 



(1 + 5Z){TTl7t2 — 7r27ri) + — (u^u — rf^d) 



(5) 
(6) 



At small temperatures, when h ^ or h = and < the original 0(4) symmetry of the meson sector breaks down 
to 0(3), which means that the a scalar field has a nonzero vacuum expectation value (vev), thus (0o) = (cr) = v ^ 0. 
Moreover, if the isospin chemical potential is sufficiently large, charged pions condensate (Bose-condensation) and 
another nonzero vacuum expectation value occurs, which additionally breaks down the symmetry to 0(2). If one 
initially introduces vev to all the four bosonic fields it can be shown that because of the remaining 0{2) symmetry 
the fields can be transformed in such a way that only two nonzero expectation value remains, from which one is the 
aforementioned v, and the other one can be either (0i) or (02), i.e. one can entirely transform out (03). Hence, we 
can choose (0i) = (tti) = p 7^ and (0i) = (-Ki) = for i = 2, 3 (see for e.g.fl^). 

Shifting the fields by their expectation values in Eq.([2]) the generating functional can be written as 



Z = / X'0X'0'X>V 



(7) 



Here 

LUn = (2?! - 



is the tree level fermion propagator matrix, which after g oing to Fourier space and introducing the 
l)/3 fermionic Matsubara frequencies, is given as (see for e.g. |l7|). 



iG: 



lMi)7o 



3 MB 



/"i)7o - liPi 



(8) 
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where pi is the spatial momentum. In ([7]) TYb = + 'Hb.2 contains all terms depending on the canonical momenta 
and bosonic fields up to second order. Other terms are contained in the £/. The non-zero condensates generate three 
point couplings, which can be directly read off from the explicit form of Ci IjA.ip . The combination 110 — TV^ can be 
written as 

n^-TYs = - 





- 6Z) 


(n^ 


-2no(i- 


hSZ)^o) 






-~5Z) 


(nl 


-203(1- 


hSZ)h) 






-SZ) 


(n? 


-2ni(i- 


h(5Z)(0i - 






- SZ) 




-202(1- 


h(5^)(<^2- 


+p))) 



(9) 

After completing into whole squares and performing the integration over the canonical momenta in one can 
identify the tree level bosonic inverse propagator matrix and from the Hnear terms the tree level equations of states 
(EoS). 

The 4x4 inverse boson propagator {G^ ) in Fourier space with the introduction of the ujn = 2?™/? bosonic 
Matsubara frequencies splits up into the inverse propagator of tt^ and a 3 x 3 coupled inverse propagator matrix of 
the remaining three bosonic fields tti, 7r2, a. In the 7ri_2 sector we switch to the charged tt^ = (tti ^ 17^2) /V2 base and 



-1 



use the following relabehng i,j = 1,2,3,4 — > 7r+, tt , ct, tts. In this way the non zero elements of G^j are 



iGl^ = i-ico,-,)^ - El^, (10) 

^(-icc;„ - Aii)2 - - Ap2 -Ap2 -V2Xvp 

iGli ' = I -Ap2 (-ic^„ + Mi)2-i?2^-Ap2 -V2Xvp 

-V2Xvp -~V2Xvp {-iuJnf - El - 2Xv^ 



■^3 

where k,l = 1,2,3, E-„^ = -h m^,^ , and the tree level tts mass square is m^^ = + A(u^ + p^). 

As mentioned earher in ([9]) the coefficients of the terms linear in the corresponding fields 9^0 (cr) and 9^1 (tti) 
determine the two non-trivial equations of states: 

EoS*;*'*' = w(m2 + A(t;2 + p2)) - = 0, (11) 
EoS*;;*' = p{m^ + A(v2 + p2) _ ^2(^ ^ ^ Q (^2) 

In the second EoS the wave function renormaHzation constant SZ is explicitly written — which unavoidably occurs 
when field renormalization is introduced followed by the canonical way of introducing chemical potentials — even if 
it is a higher order term, because we would like to emphasize the importance of its presence. From the one-loop 
fermionic contribution of the second EoS, as it will be seen later, a p\ proportional divergence arise, which has to 
be absorbed into some infinite counterterm, and this counterterm will be the SZ renormaHzation constant. In other 
words the lack of SZ would result in an uncanceled divergence in the second EoS. 

The infinite parts of the counterterms can be obtained by requiring the finiteness of the perturbative A^-point 
functions (in our case the propagator and the four point boson vertex) in the symmetric phase (f = p = 0) at 
T = pb = f^i ^ 0. Practically it is easier to obtain the infinite parts from the one-loop EoS (see Sec. 
Accordingly, the following counterterms can be found up to one-loop order with cutoff regularization. 



Sm^ = -&X{K^-m'\n^) + %N,K\ 

It in'' 



a = 12A2ln|-|L^.ln^, (13) 

where A is the three dimensional cutoff in momentum space, while lb and If are the bosonic and fermionic renor- 
malization scales, respectively. These counterterms cancel by definition all divergences at one-loop level and at 
V = p = T = pb = fJ-i = and in the broken phase the same counterterms can be used with a sHght modifica- 
tion due to the necessary resummation (see Sec. IIIip . The finite parts of Sm'^, SX and SZ are determined by the 
parameterization of the model in the broken phase (see Sec. lIVp . 
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III. EQUATIONS AT ONE LOOP LEVEL 



As it is well known from finite temperature field theory, tree level mass squares can become negative in the broken 
phase as the temperature increases. Accordingly some sort of resummation of the perturbative propagator is needed 
(see for example in [IBl)- We perform this using the optimized perturbation theory j23|, where a temperature (and 
chemical potential) dependent mass term is introduced in the Lagrangian and the difference between the original and 
the new mass parameter is treated as a (finite) higher order counterterm, 

H^ass = = -iM2(T,/^)02 _ lAm2(T,/^),^^ (14) 

Here Am^ is the finite "one-loop level" counterterm. The new mass parameter is determined by requiring that the 
inverse one-loop level tts propagator at zero external momentum stays equal to its tree level value (fastest apparent 
convergence FAC) 

M^^ = + Sm" + (A + 5\){v'^ + p^) + E^, {lu = p = 0, ,T, ^l) + Am^ = + \{v^ + p^) = ml^ , (15) 

where we indicated the resummed mass dependence of the self energy (S^g). Strictly speaking the OPT at one-loop 
order replaces m with M(T, ji) also in the internal tree level propagator lines. Since 773 is not a mixed state, M(T, /i) 
can be expressed through the tree level tts mass, which equals its one-loop level value by the condition (fT5|) . and thus 
the diagonal part of the boson propagator fTO|) can be written as a function of the tts mass instead of M{T, p), 



iGn = (-it^„ - /ii)2 - p2 - - Ap^, = (-icj„ + //i)^ - - m^^ - Ap^ 

iGl,-' = (-i^„)2 _ p2 _ ^2^^ _ 2Xv^, iGl^-' = (-ia;„)' rnl^- 



(16) 



Here and in the following the tts mass is denoted by m^j due the PMS relation (m^g = M.^^). In this way M{T,p) is 
eliminated from (fTS)) and the resummed tts mass is determined by the equation 

mi^ (T, p) = + Sm^ + (A + SX) {v^ + p^) + T.^, {uj ^ p = 0, m\ , T, /i), (17) 

where now depends on m,.^^ through the "resummed tree level" propagator matrix (flGl) . This means that the above 



procedure makes the 713 propagator (6*44) selfconsistent at p = 0, while the mixed sector of the boson propagator is 
partially resummed due to its dependence. 

The self energy in (fT7| contains bosonic as well as fermionic loop integrals. 



I],3(l^ = p = 0, m2^,T,^) = A + A^^Tr{r''G''(a;„,p,w)} 

+ ^p^^fc; (G'j'^n, P, w)G5,(t^„, p, w) + G^,(c^„, P, w)G5,(c.„, p, w)) 

+ 5F^pGf;(w„,p,//i,^B)S/jG^fc(Wn,P,W,MB)Sf;, (18) 

where T^, B^'^ denotes the coupling matrices which arise from ([9]) and are listed in the appendix. It is worth to 
note that the divergences of the self energy are cancelled by the infinite perturbative counterterms if one replaces 
TO with M(T,p) in (fT3|) . In this case, it seems that the counterterms are temperature and/or chemical potential 
dependent, however it can be proved order by order that all T and/or p dependent infinities are cancelled by higher 
order contributions (see for e.g. in jl^l and [3|). 

The scalar condensate v is determined by the vanishing of the one-loop level one point function of cr (EoSo- at 
one-loop level). 



vim^ + Sm"^ + (A + 5X){v^ + p^) + A 



^^TT{H'G''{u;n,p,pi)}+gF^Ji-{HfGf{LUn,p,pi,pB)}^ =h, (19) 

where H^'^ can be found in the appendix. Comparing (fT8|) and l|19p one can recognize a Ward identity which connects 
the symmetry breaking external field with the propagator of tt^ at zero external momentum. 



vmi^ = h. (20) 

This relation is a consequence of the remaining 0(2) symmetry, which is an axial vector rotation around the third 
isospin axis from the point of view of the chiral symmetry. Moreover (|20l) guarantees the Goldstone theorem for this 
degree of freedom (neutral pion) . 
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The pion condensate is determined through the tti one point function, which is the EoStti one-loop level, 
p (^m^ + Sm^ + (A + SX){v^ + p^) - + 5Z) + A Tr{ i?''G'' K, p, m)} 
+9F^^^r{R^G^iL0^,p,^il,flB)}^ =0, (21) 

where and R^'^ can be found in the appendix. Here the infinite part of dZ just cancels the p,i dependent divergence 
of the fermion loop integral. 

The calculation of the loop integrals in ifTS]) and (pT|) requires the diagonalization of the propagators and the 
corresponding transformation of the coupling matrices due to the non-diagonal matrix elements in (flOl) and 
The diagonalization itself is a straightforward calculation, however the eigenvalues of the boson propagator (fTO| are 
non-rational functions of Un- Thus performing the Matsubara sums is a very complicated task and it is beyond the 
scope of our paper. In order to avoid these difficulties the diagonalization was accomplished only for small p. As one 
shall see in the next section this method leads to ordinary Matsubara frequency dependence in the propagators. For 
that very reason we restricted ourselves to a Landau-type analysis of the pion condensation in a small vicinity of the 
phase boundary. Up to second order in p the mass equation (fTTl can formally be written as 



m 



2 -m2 + A«2+t(")(m2,,«,T,/ii,B) + (A + <(2)(m2^,t;,T,Aii,B))p', (22) 



T3 



while lfT9|) goes over into the simple form of (|20| due to the Ward identity. At the same order in p (pT|) can be 
rewritten as. 



pj~m^- Xv^ - r(°) (m^ , ^,, T, ^.j^) - (A + r^^) (m^ , v, T, ^^13))^' + O(p') 



(23) 



By virtue of the above equation the pion condensate may have non-zero value only if the roots of the expression in 
the square bracket are real. Assuming that A + r^^^ > 0, Eq. ((23l) yields 




— Au^ — rC') (toJ ,?j,T, utr) 



if /ij — TO^ — Aw2-r(0) > 0. This means that in this case the transition is of second order supposing that the coefficient 
of the fourth order term in l(23|) is negative. Moreover if A + r^^^ < and p?j — m? — Au^ — r^^^ < (keeping that the 
fourth order term is negative) it can be seen that the equation can have two nonzero roots, which suggests first order 
phase transition. For the calculation of the t^^'^^ and r^^'^^ coefficients in (|22|) and l(24|) diagonalization of the boson 
and fermion propagator up to order p^ is needed, which will be presented in the next section. 



IV. DIAGONALIZED PROPAGATORS FOR SMALL p 

In our approach, as was discussed previously, the next step is to determine the eigenvalues of the propagator matrices 
for small p values, in other words to find the propagating eigenmodes with help of some suitable linear transformation 
of the original fields perturbatively in p. This step is not necessary if one would like to calculate the effective potential 
(see for e.g. In case of the bosonic propagator the transformation matrix up to 0{p'^) is found to be, 

/ l-|a|V^ h{l-2av)p^ -V^ap \ 
OB=\h*{l-2a*v)p^ \-\afp^ ~^/2a* p + ©(p^), (25) 
\ V2ap V2a*p l-2|a|^pV 

where a = a{ujn,p) ~ Xv/{p? + 2\v^ + 2i(jj„/^) and b — b{ujn,p) ~ iA/(4/xw„). As it can be checked Ob is not a 
unitary transformation and it is important to note that Ob depends on the Matsubara frequency w„. With this 
transformation Ob ■ (iG^+ ^_ ^ ) • Og^ = diag(iG'~+ , iG"! , iG"^) + ©(p^), where the tilde reminds us that these 
propagators belong to the transformed (propagating) particles. It is worth to note that the new tt"*" and 7r~ particles 
are no longer charge conjugates of each other, which is a natural consequence of the presence of the pion condensate. 
After calculating the inverses perturbatively, the transformed 7r+, 7r~ and a bosonic propagators are given by 

• ^ ^ 1 2 KM + 2A^;^ - 4i^i^„) 4 

(^„ + ipif +EI P ((a;„ + i/ii)2 + ElYipj + 2\v^ - 2iwc^„) ^ ' 
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while the 713 propagator is 

In case of the fermionic inverse propagator matrix the diagonaHzation must be performed cautiously due to the 
presence of the non commuting Dirac matrices. The clearest approach is to solve the equation Oi?(iG'j^^)Oj^^ = diag 
for Of directly. In this way Of is found to be 

where fcg = (— iw„ + \i-l-b)io and the matrix is hermitian. After performing the inverse perturbatively the fermionic 
propagators are given by 

iG„/d = p —- 70-7 (31 

Pu/d - "V 8fco Pu/d - "V Pu/d - 

where = (-iwn + Mu/d)7o - 7»Pi and = ^b/3 ± fJ.i/2. 

In the appendix it is shown that all integrandus appears in the one-loop equations ifTSj) . (fT9| and ([2T1) can be 
written as traces over flavor space (in case of fermions traces also concern Dirac indices) . In this way one can insert 
the bosonic/fermionic transformation matrices given in l(25l) and (|30|) under the traces and transform the propagators 
within into diagonal form, which will lead to transformation of the corresponding coupling matrices. Thus for instance 
the trace in (|A.4p can be written as 

Tr{S'''G''} = Tr{B'''Og^OBG^Os^OB} = Tr{OBB'''Og^G''} = Tr{B^G''}, (32) 

where & = diag(G^+, G^^ , G^^), and B'^ is the transformed coupling matrix, which depends on w„ as was mentioned 
earlier. Similar expressions can be derived in case of fermions. 



V. THE PARAMETERIZATION 



Before calculating at flnite temperature and non-zero chemical potentials one has to parameterize the model at 
T = ni.B = 0. We closely follow the method presented in (2ll. [2^. [27j . Since p = at = 0, there are flve parameters, 
namely m^. A, gp, h and v, which can be fixed by setting four physical quantities — namely the pion, sigma, u and d 
quark masses plus the pion decay constant (through the PCAC relation) — to their physical values and by requiring 
the fulfillment of the l(20|) equation of state. At T = //i^b = and p = 0, the one-loop level 773 inverse propagator can 
be written as 

i(Gi;'°°P)-i = (1 - 5Z)p^ - m2 - Sm^ - (A + SX^ - J:^,{p^ = 0) ^ p^'KAP^ = 0) - S(p2), (33) 

where = dY^TrJdp^ and Sttj ~ 0{p'^). Fixing the physical 773 mass (Af^r — 138 MeV) through the 773 one— loop 
level propagator at = and using the (fT7|) mass resummation equation at p = T = pi^-Q = one obtains 

Ml = n? + Xv^ + 3XiTl;{AU, h) + T^{m„,h)) + 2glNcTl^ {ruf ,lf), (34) 

where = + 2Ai;^ and JJ^^-^'^ denotes the bosonic/fermionic tadpole integrals at zero T and /^i,b- In l(34|) 
the lack of bubble integrals is due to the fact that at p^=0 they reduce to a linear combination of tadpoles. These 
tadpoles are finite on account of the 5m? and 5X counterterms ifTS]). Note that new divergent terms do not appear in 
the tadpoles at finite temperature and/or non-zero chemical potentials, as it should be. The presence of 5Z in l(33|) 
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FIG. 1: Bosonic renormalization scale dependence of the fermion renormalization scale If, the mass parameter (left panel) 
and of the coupling A, the tree level a mass rricr (right panel). 



is also important because this term renders the dependent part of the propagator finite. It involves also a finite 
renormalisation: 



(p2 = 0) 



Ml - mt 



.log($) 



167r2 



{ml - Mif 



167r2 



log 



(35) 



(36) 



and the one-loop level PCAC relation depends on &Z^'^ as follows 

«[iGif°P(p' = 0)]-i=/.A/2(l 

where /^r = 93MeV is the pion decay constant. Since 5Z^^ depends on the fermionic renormahzation scale its actual 
value is tunable. Thus we required the vanishing of 6Z^^^ to fix the fermionic renormalization scale. This requirement 
makes the PCAC relation simpler and thus the values of v, gp, h offer themselves immediately 



/it J 



gp = 2^, h = UMl 



(37) 



where m/ = 938/3 MeV is the constituent u,d quark mass. At this point v, gp, and h are known and can be 
expressed from l(34|). 

The remaining unknown A parameter is determined by fixing the physical a mass at one-loop level and at zero 
external momentum, that is 



Mi 



5m 4 



3(A + a)i>2 +S^(p2 = 0) 

+ ?,\{T^{M^,h)+T^{m,,k)) + l%X^v''Bl{m,,k) + &\^v''Bl{M^,h)+QglTf^''{mf,lf), (38) 



where Bq and t/''^ come from the temperature and chemical potential independent part of the bosonic/fermionic 
bubble diagram at zero external momentum and degenerate masses. It is worth to note that the same infinite 
counterterms render the above equation as well as the equation of M^r finite. 

We use Mfj = 500 MeV for the physical <t mass. This choice seems somewhat arbitrary, because the a meson is 
a broad resonance rather than a particle with well-defined mass (see [l^l and references therein). Thus it would 
be more appropriate to identify the mass and width of the a meson through the pole of its spectral function [29| . 
However, we checked that varying M^ in the 400 MeV - 750 MeV range produces just the same order of uncertainties 
in the thermodynamical results as the variation of the lb boson renormalization scale itself. Accordingly, it is enough 
to analyze the h dependence of the different parameters, hence here and in the following we present our results which 
correspond to the choice M^ — 500 MeV at zero T, ^lb. On the left panel of Fig. [T] the bosonic renormahzation scale 
dependence of m^ and while on the right panel the tree level cr mass and the A parameter are shown. As one 
can see on the right panel, the tree level a mass equals its one -loop level value at Ih ~ 600 MeV thus the a mass 
becomes selfconsistent at this point (but only for zero T, /Ui^bO- addition m and If moderately depend on the 
renormalization scale around this point (see left panel of Fig. [l|) . According to the above arguments we choose the 
following scale range: lb € [400 MeV, 800 MeV] for the thermodynamical calculations. 
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FIG. 2: Temperature and chemical potential dependence of the scalar condensate v. Left panel shows v{T) at Ati,B = and 
h = 400, 600, 800 MeV. Right panel shows vi^is) at m ^ and T = 20, 50, 80, 100 MeV. 



VI. RESULTS AT LOWEST ORDER IN p 

As was derived in Sec. IIIII the expectation value p of the pion field is determined by l(24|) up to second order in 
p. From this equation it can be immediately seen that the second order boundary for the occurrence of the pion 
condensation in the pi — pB — T space is determined by 

pj-m^~ Xv^ ~ r(") (rr^ , v, T, ^ij^) = 0. (39) 

Moreover, the pi^B and T dependence of v and are determined by l(20| and l[22l) . which have to be solved at 
p = 0. At p = (|22)) has the same form as l(34|) with the slight difference that now the tadpoles which appear therein 
have to be calculated at finite temperature and chemical potentials, while the form of (|20| is unchanged. Using the 
explicit expression of r^°^ , Eq. ((39l) can be expressed as 

pj ~ ml^ (T, Ml , MB ) - (T, MI , MB ) = 0, (40) 

where B}-^°°p is the remaining part of r^°) after subtracting from it mj^ . It contains one-loop bosonic and fermionic 
contributions. From l(40|) it is obvious that at one-loop level the condensation does not start exactly at p\ = rriT^^, 
as it is commonly expected, but it is shifted to some extent by i?^"'°°P. It's worth to note that the deviation iji-'oop 
does not vanish identically even if the resummed pion mass is defined as the pole of the propagator. 

First, we investigated the temperature and chemical potential dependence of the scalar condensate v at mi = 
and different values of h by solving l(20|) and (|22l) , which can be seen in Fig. [21 On the left panel the temperature 
dependence of v is shown at h = 400, 600, 800 MeV and at zero chemical potentials. As it is expected, the scalar 
condensate shows a smooth crossover as the chiral symmetry is restored at a pseudocritical temperature (Tc) around 
150 MeV, which is in good agreement with the continuum limit lattice result found in [sOl- Moreover, the v{T) curve 
and consequently Tc slightly depend on 4- Hence, in the forthcoming we use the fixed scale h = 600 MeV. On the right 
panel the baryochemical potential dependence of v can be seen at 4 — 600 MeV scale, and at T = 20, 50, 80, 100 MeV 
temperatures. At small temperatures the transition is of first order, while for large temperatures it is of analytic 
crossover type, indicating the existence of a critical endpoint (CEP), where the transition changes from first order 
to crossover with increasing temperature (see e.g. in [2ll] and references therein). As one can see on Fig. [2] the 
temperature at the CEP is around 50 MeV, which is much lower than the lattice result presented in Ref. j3l|, however 
this is a common feature of effective models (for two/three fiavors see e.g. Refs. [2^, [2l|). The critical/pseudocritical 
baryochemical potential values range from ~ 600 MeV to ~ 1000 MeV depending on the temperature. 

Next, solving ((20l) and ([22]) for different T, mi,Mb values and tracking the fulfillment of the l(40|) condition we 
determined the second order critical surface of pion condensation, which can be seen on Fig. [31 As a function of the 
isospin chemical potential at fixed mb the region of pion condensation starts very steeply at around 130 MeV. After 
that steep jump a plateau starts, which decreases slowly towards higher values of mi- Moreover, it can be seen that 
the region of pion condensation shrinks with increasing mb and it even disappears at around mb = 830 MeV (and 
Mi = 131 MeV), a behavior which is in accordance with previous effective model studies (see e.g. Sec. VILA of [l3|)- 
This is understandable physically, since at large mb the condensate p is basically determined by the difference of the 
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FIG. 4: Phase boundary of pion condensation and the /ii = condition as function of the isospin chemical potential for 
different /ib values. 

u — d quark contributions in the EoS, and this difference is decreasing with increasing fiB at fixed /xi, because the 
Fermi-Dirac factor in the contribution of u depends on + Mi j while in the contribution of c? on /xb — Mi ■ It is worth 
to observe that on the surface a gradually increasing missing part starts from about /ib = 415 MeV and /ii = 221 MeV. 
In that region the M/ ~ ™^ ~ -^^^ ~ ^^^^ combination (numerator of Eq. (|24| ) is negative, which means that the 
transition is not second order anymore. Strictly speaking if in this region there is a nonzero solution to p this can 
happen only if A + < (denominator of Eq. l(24|)) and in this case the transition is of first order. 

To analyze the surface in detail, two sections taken at mb = MeV and mb = 400 MeV are plotted in Fig. [4]together 
with the Mi = ™7r3 curves. At both baryochemical potentials the condensation starts at around mi = 131 MeV, which 
is slightly below the m^., = 138 MeV pion mass, this deviation is due to the corrections iji^'oop in ||40p . Moreover, 
the deviation widens as mi increases. Another interesting feature is that at fixed high temperature as we increase the 
isospin chemical potential the condensation evaporates above a certain m value, which is in accordance with Ref. [lB |. 
where this phenomenon was observed in case of two flavor PNJL model. In the region where the condensate already 
evaporated (p = 0) the chiral symmetry is almost totally restored {v ~ 0). 

Finally we calculated the one-loop pole masses of the charged pions on different temperatures as a function of the 
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FIG. 5: /ii dependence of the one-loop pole masses of charged pions at different temperatures (T = 90, 140, 150 MeV). The thin 
lines indicate the condensed phase for each temperature (where our calculation is not valid). 

isospin chemical potential at p = 0. For this we calculated the self energies of the charged pions, then we solved the 
following pole mass equations: 



where we have implicitly used the already known solutions f(r, /ij^) and TOttj (T, /.tj b). The one-loop masses are 
plotted on Fig. O On the figure thick Hues represent the sectors, where our calculation is vaHd, that is p = 0. At the 
thin line parts of the curves the pions are condensed, thus we should go beyond the lowest approximation in p to get 
correct results in that sectors. It is worth to note that when the condensation sets in none of the one-loop charged 
pion masses becomes zero as opposed to the masses defined through the dispersion relation (see e.g. Another 
interesting thing is that at larger values of where the condensation has already evaporated and our calculation is 
valid (thick line parts) the charged pion masses are still different, however p = 0. This difference is due to the fact 
that the isospin background acts differently on the charged pions (with opposite signs). 



In this paper we studied the pion condensation in the framework of SU{2)l x SU(2)ii constituent quark model 
with explicitly broken symmetry term in the presence of baryochemical potential. The model was parameterized 
at one-loop level and optimized perturbation theory was used for the resummation of the perturbative series. The 
one-loop equations were expanded in powers of p, and a Landau type analysis was performed for the phase boundary 
at lowest order in p. A simple condition for the boundary of the pion condensation was set up, and we argued that 
this condition gives a second order surface in the /zi — — T space. The temperature and renormaHzation scale 
dependence of the scalar condensate v was investigated, and a mild renormaHzation scale dependence was found. 
At zero baryochemical potential the pseudocritical temperature is in accordance with results found on lattice [30| . 
Using the condition for the pion condensation the second order surface was determined. It was found that the surface 
starts steeply with increasing pi at fixed pB and towards large values of pB the pion condensed region shrinks and 
even disappears at around ^b = 830 MeV. However, at such a high energy one should take into account the effects 
of the strange quark. Investigating different sections of the surface it was showed that at one-loop level the pion 
condensation curve sHghtly differ from the pi = m^j curve at small pi and this deviation increases with increasing 
Pi. We also studied the dependence of the charged one-loop pion masses on the isospin chemical potential. As a 
continuation of the present study the analysis can be extended to higher order in p, with which for instance the scaling 
properties around the surface, the dispersion relation at one-loop level and diflerent phases of the condensed matter 
(BCS/LOFF) could be investigated. 




(41) 
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APPENDIX: COUPLINGS 



The coupHng matrices appearing in ifTSj) . ifTO]) and (|2l|) can be obtained from the interaction term of the shifted 
fields in the Lagrangian 

^'^'^ + ^'4'Ta4>a'4> > ^ ((tTi + p)^ + + nl + {(J + v)'^)'^ + (To(cr + V) + i75(Tl(7ri + p) + T2-K2 + TsTTs)) ?/). 

(A.l) 

The four and three point couplings of ijA.ip determine the coefficients of the tadpole and bubble terms of (fT8| . 
Including the symmetry factors of the corresponding graphs these are: 



/I o^ 

10 

10 

VO 3/ 



/O p \ 

p 

V2v 

\p p ^/2v J 



-75 
75/' 



(A.2) 



where in case of bosons the same convention is used for the labeling of the matrix elements as in (fTO|) . The bosonic 
bubble contribution in (flSl) can be rewritten as 



J2 B\^BUg\iG], + G]iG\,) = 2Tr{G''^B^G''B'}, 



(A.3) 



and by virtue of the structure of B^ and G'' one can factor out the tt^ propagator from the above expression as follows 



TrjC"' B'G^B"} = G^ TrjS" G"}, 



(A.4) 



where 



B" 



( p^ V2vp 0\ 

V2vp V2vp 

V 0/ 



The fermionic bubble contribution can be formulated as, 



Y: GiB(P^,Bi = Tr{Gf^BfGfBf] = Tr{BfGfBfGf^l 

i,j,k,l 



(A.5) 



(A.6) 



where the trace is over fiavor as well as Dirac indices. 

Moreover, in equations of states lfT9|) and (|2T|) the tadpole coefficients are also determined by the three point 
coupHngs of (|A.1[) . 



( 2p p v/V2 Q\ 
p 2p 
^l/^/2 p 

V Op/ 



2 lo 1 



(A.7) 



and 



I V p/V2 Q\ 

V pI^/2 

p/^/2 p/V2 iv 

V Of/ 



2 1 75 



(A.8) 
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